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ABSTRACT 

We use four Chandra gratings spectra of the neutron star low-mass X-ray binary 4U 1820-30 to better 
understand the nature of certain X-ray absorption lines in X-ray binaries, including the Ne II, Ne III, Ne IX, 
O VII, and O VIII lines. The equivalent widths of the lines are generally consistent between the observations, 
as expected if these lines originate in the hot interstellar medium. No evidence was found that the lines were 
blueshifted, again supporting the interstellar medium origin, though this may be due to poor statistics. There 
is apparent variability in the O VIII Lya line equivalent width providing some evidence that at least some of 
the O VIII absorption arises within the system. However, the significance is marginal (2.4-a), and the lack of 
variation in the other lines casts some doubt on the reality of the variability. From calculating the equivalent 
hydrogen column densities for a range of Doppler parameters, we find they are consistent with the interstellar 
origin of the lines. Additionally we fit the spectra with photoionization models for locally absorbing material, 
and find that they can reproduce the spectrum well, but only when there is an extremely low filling factor. 
We conclude that both the ISM and local absorption remain possible for the origin of the lines, but that more 
sensitive observations are needed to search for low-level variability. 

Subject headings: accretion, accretion disks — stars: neutron — X-rays: individual (4U 1820-30) — X-rays: 
binaries — X-rays: ISM 



1. INTRODUCTION 

High-resolution X-ray spectroscopy from Chandra and 
XMM-Newton allows for the study of absorption and emis- 
sion lines in Galactic X-ray binaries. While emission lines 
are clearly associated with the source s themselves (e.g., 
ICottam et all 1200 lalbt ISchulz etai1l2001l) . that is not always 
the case for absorption lines. In some systems, the absorp- 
tion lines are also clearly associated with the sourc e, as they 
are observed to be variable and/or blueshifted (e.gjLee et all 



20021: iParmar et alj2002l;IUeda et alJ2004tlBoirin et alj 2005t 
Diaz Trigo et all 120061: IMiller et al.M2006allbh . However, in 
other sources the lines are interpreted as being associated 
with absorption by the hot ionized interstellar medium, which 
is thought to be present in the disk an d halo of the galaxy , 
with temperatures around 10 6 K (e .g., [Fu tamo to et alJl2004t 
lYao & Wang|l2005llJuett et al.ll2006l and references therein). 

From studies of the absorption lines in the black hole X- 
ray binary GRO J 1655-40, the plethora of highly-ionized 
lines that are detected are all seen to be blueshifted - ev- 
idence for a mass outflow into our line of sight, presum- 
ably from the accretion disk. In fact, detailed modeling of 
the lines in this system sugge st that magnetic fie lds may 
be the force driving the wind dMiller et al.l l2006ah . Other 
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black hole systems also show evidence for highly-ionized out- 
flows such as H 1743-322 GX 339-4 an d GRS 1915+105 
(IMiller et al.ll2006bl 12004 iLee et alj[2002b . Disk winds are 
also observed from the accretion disks arou nd supermassive 
black holes in Active Galactic Nuclei (e.g . Crensha w et all 
Il999l iKaastra et al.ll2000l iKaspi et alJl2000h . 

Several neutron star low-mass X-ray binaries also display 
clear evidence for absorption local to the source. For ex- 
ample, GX 13+1 shows blueshifted absorption li nes indicat- 
ing a n outflow with a velocity of ^-400 km s" 1 (Ue da et ail 
2004). Moreover, in the X-ray dipping sources, strong absorp- 
tion lines are observed as these objects are inclined close to 
the orbital plane (e.g.. ISidoli et alJl200lt: IParmar et alJl2002b 
iBoirin et alj2005l:IChurch et alj2005l:lDfaz Trigo et al.l2006l) . 
On the other hand, the nature of the absorption lines ob- 
served in many other neutron star X-ray binaries is less clear 
- the lines are weaker and there is no evidence of signifi- 
cant blueshifts, supporting the interpreta tion that these lines 
are due to the in terstellar medium (e.g.. lYao & Wangll2005l : 
Uuett et al]l2006h . 

As part of the Chandra HETGS Z/ Atoll Spectroscopic Sur- 
vey (CHAZSS) we have obtained sensitive observations of 6 
neutron star low-mass X-ray binaries at 2 separate epochs us- 
ing the Chandra high energy transmission grating spectrom- 
eter (HETGS). The multiple epochs allow for a search for 
variability in the absorption lines present in these systems; 
this is key for learning about the nature of the absorption 
lines. In this paper, we present the observations of the sys- 
tem 4U 1820-30 from this survey. Tight limits can be set on 
the size of any absorbing region associated with a binary if 
the system is ultra-compact in nature. Therefore, the ultra- 
comp act nature of 4U 1820-30 dStella. Priedhorskv, & White! 
119871) makes it a particularly interesting source to study ab- 
sorption lines. 
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2. PREVIOUS CHANDRA OBSERVATIONS OF 

4U 1820-30 

4U 1820-30 is an accreting neutron star low mass X- 
ray binary in the globular cluster NG C 6624, which has a 
well determined distance (7.6 ± 0.4 kpc Kuulkers et al. 2003) 
and reddening (E( B-V) = 0.32 ± 0.03 Kuu lkers et a7Tl2003t 
iBohlin et al.lll978l) . This reddening is consistent with the to- 
tal HI column density of our Galaxy in this direction deter- 
mined from 21 cm radio o bservations, Nh = 1.5 x 10 21 cm" 2 
dDickev & Lockmanll 19901) . We adopt the above column den- 
sity throughout this paper. 

Absorption lines and edges in the X-ra y spectrum of 
this so u rce have been studie d pr eviously by |Futamoto et all 

(2004) , lYao & Wangl (120051) and lYao & Wa n gl (120061) using 
an LETG/HRC Chandra observation (Obsld 98), and by 
Juett et al.ld2004l).lYao & Wang (20051). lJuett et al.ld2006l) and 
Yao & Wangl (120061) using HETG/ACIS Chandra observa- 
tions (Obslds 1021 and 1022). A summary of the Chandra 
observations of this source is given in Table [TJ We briefly 
s ummarize the main resu lts of these previous studies below. 

iFutamoto etail (120041) clearly detect the O VII Hea, O VIII 
Lya, and Ne IX Hea absorption lines in the LETG spectrum, 
and put upper limits on the O VII He/3 line. From their curve 
of growth analysis and photoionization modeling, they deduce 
that all oxygen will be fully photoionized if the absorbing col- 
umn is located close to the binary system, and therefore at- 
tribute these lines to hot gas in the ISM. However, we note 
that they make two tacit assumptions - that the filling factor 
of the gas is 1, and that all the absorption arises from the same 
gas. 

The studies by Uuett etail (120041) and Uuett et all (120061) are 
concerned with the absorption edges from oxygen, neon and 
iron, comparing t he spectra of a number of X-ray binaries. In 
Uuett et all d2004l) the structure around the oxygen edge in 7 
X-ray binaries is studied, looking at the O I, O II, and O III 
absorption lines in that region. The edge and thes e low ion- 
ization lines are attributable to ISM absorption. In lJuett et al.l 
(200$ the neon and iron L edges are studied, and the nature 
of the neon absorption lines discussed. These authors come 
to the conclusion that the neon lines in the low-mass X-ray 
binaries studied are consistent with predictions for the ISM, 
except for the N e IX line detected in the black hole X-ray bi- 
nary GX 339-4 iMiller et al] (120041) which shows a column 
density significantly higher than predicted for the ISM, con- 
firming that the largest contribution to this line (in this source, 
at least) is from local absorption. 

From a study of 10 low-mass X-ray binaries. lYao & Wand 

(2005) find that the detected Ne IX, O VII and O VIII absorp- 
tion lin es are consistent with the hot ISM origin. lYao & Wangl 

(2006) focuses specifically on 4U 1820-30, discussing the 
implications of the Ne and O abundance in the ISM from the 
absorption lines/edges in the Chandra observations. 

3. OBSERVATIONS AND DATA REDUCTION 

All the previous observations discussed above suffer from 
pileup, where a high count rate leads to more than one photon 
arriving at a pixel per read out. This affects the shape of the 
spectrum, typically making it harder as the sum of the ener- 
gies of the multiple photons that arrived between each read- 
out is detected as one count. It is unclear how pileup affects 
the measurement of absorption lines. It likely modifies the 
continuum level and shape, which would in turn change any 
measured equivalent widths (which depends on the continuum 



flux across the line). Currently, there is no unambiguous way 
to model pileup in Chandra gratings spectra. To minimize 
pileup we operated Chandra with the High Energy Trans- 
mission Grating Spectrometer in continuous clocking mode 
to significantly increase the temporal resolution. Continuous 
clocking mode allows milli-second timing at the expense of 
one dimension of spatial information. 

As part of CHAZSS, 4U 1820-30 was observed twice with 
Chandra for approximately 25 ksec each time, first on 2006 
Aug 12 (Obsld 6633), and later on 2006 Oct 20 (Obsld 
6634). The HETGS -dispersed spectrum was read-out with 
the ACIS-S array operating in continuous-clocking mode. We 
reduced these data using CIAO version 3.3.0.1 and follow- 
ing the standard analysis threads 8 . Using the tgextract tool 
spectra the + 1 and - 1 order were extracted at the nominal in- 
strument resolution for both the high-energy grating (HEG) 
and the medium-energy grating (MEG). Using the standard 
redistribution matrix files (RMFs) from the CALDB, the an- 
cillary response files (ARFs) were generated using the fullgarf 
script. We carefully inspected the MEG +1 and -1 orders to 
confirm there was no wavelength shift between the two or- 
ders. The spectra and ARFs were then added together us- 
ing the add_grating_spectra, giving a first-order spectrum for 
the MEG. Similarly, we created a first-order spectrum for the 
HEG. 

In addition to analyzing these two new Chandra observa- 
tions, we also re-analyze the two previous Chandra HETGS 
observations of 4U 1820-30 (Obsld 1021 & 1022) to allow 
for a search for variability in the absorption lines. These two 
datasets were not operated in continuous-clocking mode, and 
suffer from pile-up. We reduce the data following the stan- 
dard analysis threads, extracting the spectra and ARFs using 
the same method described above. For consistency, we only 
analysis the HETGS observations here and not the LETGS 
observation which has a lower spectral resolution. 

4. ANALYSIS AND RESULTS 
4. 1 . Absorption lines 

When searching for absorption lines, we use only the MEG 
spectra, as the MEG has a higher effective area over the re- 
gion of interest compared to the HEG. We inspected the spec- 
tra over 2A segments from 2A to 24A using the ISIS spectral 
fitting package. In each segment, we fitted a simple power- 
law (modified, where appropriate, by an absorption edge due 
to neutral elements in the ISM). Within each segment, we 
searched for any significant absorption lines present. We fit- 
ted a Gaussian to the detected absorption lines whose central 
wavelength, width, and normalization were all allowed to be 
free parameters. The equivalent width of the line was deter- 
min ed from the fitted Gaussian. F or line identifications we 
use IVerner. Verner. & Ferland d 19961) as a reference. How- 
ev er, for Ne II and N e III we use the wavelengths determined 
by lJuett et afl {2006). They find the weighted mean of the 
central wavelengths of these lines measured in multiple X- 
ray binaries. The wavelengths for these lines from different 
sources are seen to consistent with each other, and are all sys- 
tematically shifted by a small amount (~20mA) fr om the the- 
oretical value determined by Behar & Netze n d2002l) . though 
within their stated errors. The lJuett et all (120061) values are 
known most accurately, thus we use them throughout the rest 
of the paper. 

8 http :/, http ://cxc .harvard . edu/ciao/threads 
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In the two new Chandra observations, we detect several 
X-ray absorption lines. We define a detected line as one 
for which fits with a Gaussian give an equivalent width that 
excludes zero at more than the 2a level of confidence (note 
that we are being conservative by using the significance with 
which we can measure the equivalent width as the significance 
of the detection). In increasing wavelength (theoretical/best 
observational wavelengths are given in brackets), the detected 
lines are: Ne IX (13.447 A), Ne III (14.508 A), Ne II (14.608 
A), O VII (18.629 A), O VIII (18.967 A). All these absorption 
lines are detected in both observations, however, some were 
not detected in the previous Chandra observations, likely due 
to the lower sensitivity of those shorter observations. Whilst 
the O I absorption line at ~23.5A was detected in all obser- 
vations, we have not attempted to fit this line due to the com- 
plicated structure around this region due to multiple oxygen 
fines. The stru cture in th i s reg ion has been studied in detail 
previously by lJuett et al.l (120041) . The interstellar absorption 
at this wavelength also significantly increases the noise in the 
data, therefore the LETG observation of 4U 1820-30 previ- 
ously studied present the best available dataset for analysis of 
this region. We also note that the O VII line at —2 1 .6 A was not 
detected in any of these observations, due to the low signal-to- 
noise in that region. Again, the previous LETG observation 
of 4U 1820-30 gives the best detection of this line. 

In cases where an absorption line had been detected in one 
of the other observations, but not in the observation we were 
studying, a 95% confidence level upper limit for the line flux 
and equivalent width was determined. This was achieved by 
fixing the line wavelength at the theoretical value, fixing the 
width at a value smaller than the instrument resolution (as 
would be the case if the line was unresolved), and then in- 
creasing the line flux until the \ 2 value for the fit increased 
by 2.7. We justify fixing the width smaller than the instru- 
ment resolution as for the MEG, the resolution corresponds to 
a FWHM of -350 km s" 1 at 18. 967 A, and higher at shorter 
wavelengths. This is significantly higher than the expected 
thermal broadening in the ISM. Table|2]gives the line proper- 
ties determined from the spectral fitting. In Fig. Q]we show 
the absorption lines detected in the two new Chandra obser- 
vations. 

As Ne IX is detected at 13. 447 A one should also expect to 
find additional Ne IX lines at 11. 547 A and 11.000A. How- 
ever, we do not clearly detect these lines in any of the obser- 
vations. If the line is not saturated, the ratio of the equiva- 
lent widths should follow the ratio of the oscillator strengths 
times the ratio of the wavelengths squared. The mean value of 
the equivalent width of the Ne IX (13. 447 A) line is approxi- 
mately 5mA. From the oscillator strengths, this would put the 
expected equivalent widths for the Ne IX lines at 1 1 .547 A and 
11.000A to be 0.8mA and 0.3mA, respectively. In Ob- 
sld 6633, we note that there may be a weak absorption line 
present at A = 11.550 ±0.007 A , with an equivalent width 
of EW = 1 .8 ± 0.9 mA, though this is only a 2a detection. 
This measured equivalent width is consistent with what is pre- 
dicted from the oscillator strengths. The line at 11.000A is 
not detected, and we determine the 95% confidence level up- 
per limit to be 1.7 mA - the line would be too weak to de- 
tect given the noise in the data. In Obsld 6634 these Ne IX 
lines are not detected, and we find the upper limits on the 
equivalent width for the lines at 11. 547 A and 11.000A to 
be 1.6mA and 1.3mA respectively, consistent with what 



would be expected from the oscillator strengths. We do not 
determine upper limits for these lines from the Obsld 1021 
and 1022 datasets, as they are substantially less sensitive, 
and therefore the lines could not be detected at the predicted 
equivalent widths. While the upper limits do not rule out that 
the line is saturated, we cannot be conclusive. 

Similarly, we try to determine whether the O VIII lines are 
saturated by determining limits on the O VIII Ly/3 line at 
16.006A. Fitting a Gaussian, with central wavelength fixed 
at 16.006A to Obsld 6633 and 6634 we measure equivalent 
widths of 3.0 ± 1 .5 mA, and 3.4 ±2.8 mA respectively. This 
gives the ratio of the equivalent widths of the O VIII Ly/3 to 
Lya as 0.47 ± 0.32 and 0. 1 8 ± 0. 16 for Obsld 6633 and 6634 
respectively. For comparison, if the lines are not saturated we 
would expect this ratio to be 0. 14 given the oscillator strengths 
and wavelengths. The measured value from Obsld 6634 is 
clearly consistent with this ratio, and in Obsld 6633 the ratio 
is not significantly different. Therefore the O VIII lines are not 
highly saturated. 

Another check on saturation can be performed using the 
measured equivalent width of the O VII Lya line (measured 
in pre vious Chandra LETG observations by Futamo to et al.l 
2004) and comparing it with the weighted average of the 
O VII Ly/3 equivalent widths measured here. We find that 
the ratio of O VII Lya equivalent width to Ly/3 equivalent 
width is 5.5 ±2.2. F rom the line wave lengths and oscillator 
strengths (taken from lVerner et al.ir9 96). the expected ratio in 
the non-saturated regime is 6.4, consistent the ratio from the 
observed equivalent widths, suggesting that the O VII lines are 
not highly saturated. 

4.2. Column Densities 

From the measured equivalent widths of the lines, we de- 
termine the column densities. Firstly we assume that the 
lines are not saturated and on the linear part of the curve 
of growth. Under this assumption, the column density and 
equivalent width are associated via W\ = (ne 2 / M e c 2 )Nj\ 2 fij = 
8.85 x \0~ n Nj\ 2 fij, where W\ is the equivalent width (in 
cm), Nj is the column density of a given species, A is the line 
wave length (in cm) and is the oscillator strength (Srjitzer 
[19781) . 

However, as we cannot conclusively determine whether the 
lines are saturated, we also calculate the column density for a 
range of Doppler parameters, b. For a lower limit for the line 
broadening, we estimate the FWHM of the Galactic rotation 
in the direction of the source from a Besancon model of the 
Galaxy to be 50 km s" 1 . We also evaluate the column densities 
assuming b = 100 km s" 1 and b = 200 km s" 1 . The latter being 
an upper limit on the wid th of O VI emission lines from UV 
observations with FUSE (Otte & Dixonll2006l) . Given the os- 
cillator strength, line wavelength, observed equivalent width 
and assumed Doppler parameter, we determine the required 
column density for each line in Obsld s 6633 and 6 634, us- 
ing the curve of growth equations in Spitzer (1978). For the 
blended O VIII fine structure components we compute the col- 
umn densities numerically. Table [3] gives the calculated col- 
umn densities. 

In order to test whether these column densities are consis- 
tent with the interstellar medium, we converted the column 
density for each species into an equivalent hydrogen column 
density. To do this we require the ionic fraction for each ion, 
which for collisional ionization (the mechanism present in the 
ISM) depends on the temperature. The low ionization lines 
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(Ne II and Ne III) will be in a lower temperature gas than the 
higher ionization lines (Ne IX, O VII, O VIII), thus we esti- 
mated the temperatures separately. Comparing the observed 
ratio of Ne II/Ne III with the ratio of the ionic fractions for 
those ions at a rang e of temperatures (as calculated for colli- 
sional ionization by Mazz otta et al.llT998h . we determined the 
best-fitting temperature of the gas. Given that temperature, 
the i onic fractions at that temperature and the ISM abundances 
from Wilms et al.1 £2000) we converted the ionic column den- 
sities into equivalent hydrogen column densities (see TableO. 
We found the temperatures for the Ne II and Ne III gas to be 
in the range (4.8-5.2) x 10 4 K. 

Similarly for Ne IX, O VII, O VIII we determined a tem- 
perature from the r atios of O VII/O VI II and O VII/Ne IX (as- 
suming Ne/O from Wilms et al. 2000) compared to the ionic 
fractions from lMazzotta et aTl (119981) . and then used the esti- 
mated ionic fractions and abundances to convert to equivalent 
hydrogen column densities (see Table |5J. The temperatures 
we determine for this gas is in the range (1.3-2.0) x 10 6 K, 
consistent with the hot phase of the ISM. 

The equivalent hydrogen column densities that we estimate 
are all consistent with having a lower column density than 
the column to the source as seen in H I (1.5 x 10 21 cm -1 , 
iDickev & Lockmanll990l) . 

4.3. Photoionization modeling 

In order to assess the possibility of the observed absorp- 
tion lines originating from local absorption in the X-ray bi- 
nary system (rather tha n the ISM) we used the p hotoioniza- 
tion code XSTAR v2.1 dKallman & Bautistall200l to model 
the observed spectrum in ObsID 6633 and 6634. For the in- 
put ionizing spectrum we use the unabsorbed source spectra 
for each observation, scaled to the observed luminosities (see 
section l4~4l for details on determining the source flux, we as- 
sume d = 7.6 kpc giving 0.5-10 keV luminosities of 5.7 x 10 37 
and 8. 1 x 10 37 erg s" 1 for ObsID 6633 and 6634, respectively). 
For each observation we then created a grid of absorption line 
models allowing the column density and ionization parame- 
ter, £ = L/nR 2 , to vary. In our models we assume solar abun- 
dances, a covering fraction of 0.5, and turbulent velocities of 
200 km s" 1 . From these grids of models we then fit the ob- 
served spectra (again using ISIS) from 13A to 20A with a 
power-law (with Galactic absorption fixed at 1.5 x 10 21 cm" 1 ) 
convolved with the absorption line models calculated from 
XSTAR. This allows us to determine under what conditions 
the observed spectrum can be reproduced by photoionization 
of local absorbing gas. 

In Fig. |2] we show the best-fitting models from the pho- 
toionization modeling. We find that for both observations, the 
observed Ne IX, O VII, O VIII lines can be reproduced well 
with xl = 113 for ObsID 6633, and xl = 115 for ObsID 
6634. Note that the Ne II and Ne III lines are not included 
in the code, and thus are not reproduced. 

For ObsID 6633 we find the best fit with N H = (6.1 ± 1.5) x 
10 19 cm" 2 and log£ = 1.36 ± 0.05 (where £ = L/nR 1 ). For 
ObsID 6634 we find the best fit with N H = (7.6 ± 1.3) x 10 19 
cm" 2 and log £ = 1 .34 ± 0.03. Gi ven that the binary separation 
for 4U 1820-30 is - 10 10 cm (IStella et alJll987l) . with the 
ionization parameters we find the gas density would have to 
be n = 2.5 x 10 16 cm" 3 for ObsID 6633 and n = 4 x 10 16 cm" 3 
for ObsID 6634 for the largest plausible radius for absorbing 
gas within the system. The implies that the filling factor, /, 
is extremely low with / = Nn/nR = 2.4 x 10" 7 for 6633 and 



/ = 1 .9 x 10~ 7 for 6634 . Since / scales with R, smaller radii 
within the system make the problem worse. 

It seems unlikely that the gas would be in a shell, as the 
filling factor would mean that the thickness of the shell would 
be very small (on the order of only tens of metres for gas 
within the system). Additionally it is also unlikely to be a 
shell covering the entire 4tt or else one would also observe 
emission lines. One possibility is that thermal instabilities 
produce small, dense structures, or alternatively there could 
be partial covering by larger blobs, which could be driven by 
a wind. Therefore, while it seems possible that locally pho- 
toionized material could account for the abso rption, it does re- 
quire f airly extreme parameters. We note that Futamoto et al. 
(2004) use CLOUDY to run a numerical simulation for pre- 
vious observations of 4U 1820-30. These authors concluded 
that photoionized gas within the binary system could not re- 
produce the observed line column densities. In summary, the 
ISM may be the more likely source of the absorption lines, but 
our models demonstrate that photoionization cannot be ruled 
out, especially given the presence of variability. 

4.4. Source flux 

In order to compare absorption line properties with the 
source properties, we determine the source flux from the best 
fitting continuum model to the HEG s pectrum. For th is broad- 
band spectral fitting we use XSPEC dArnaudlll996l) . We do 
not fit the MEG and HEG simultaneously due to the signifi- 
cant pileup in the first two observations, where the observa- 
tions were not operated in continuous clocking mode. As the 
MEG has a larger effective area than the HEG at low ener- 
gies (by about a factor of 2) it is more substantially affected 
by pileup. This is clearly apparent when trying to fit a con- 
tinuum model to just the MEG spectra in these observations. 
There is an excess of counts at high energies and deficit at 
low energies, as expected by pileup. The HEG is much less 
affected, and thus we choose to fit to just the HEG 1st order 
spectra. For consistency, we use this method for fitting to the 
new Chandra observations also, even though they do not suf- 
fer from pileup. 

We choose to fit an absorbed blackbody plus power-law 
model to the HEG spectra over the energy range 1.2-8 keV 
with the column density is fixed at the galactic value towards 
this globular cluster for all observations. Model s to neu- 
tron s tar X-ray binary spectra can be degenerate dLin et al.1 
120071) . especially over a short energy range, such as covered 
by the HEG. The model we adopt is fiducial and reflects a 
physically-motivated scenario consisting of thermal and non- 
thermal emission. It allows us to characterize the source flux 
well because it provides a good fit to the spectrum. Impor- 
tantly, it is a simple continuum model and so is easily repro- 
ducible, but, due to modeling degeneracies the specific pa- 
rameter values derived should be regarded with caution. The 
fitted spectral parameters are given in Table |4] 

We determine the 0.5 - 10 keV source flux, by extending 
a dummy response to lower and higher energies. We deter- 
mine the uncertainty in the flux by propagating the errors on 
the individual parameters. The calculated fluxes are given in 
Table [4] We note that the uncertainties in the fluxes may be 
dominated by calibration uncertainties, and could be closer to 
around 6% (the difference in flux we find between fitting to 
the MEG and HEG) than those we determine from spectral 
fitting. 

We check whether pileup alters the determined fluxes sub- 
stantially by comparing the calculated source fluxes with the 
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1-day average RXTE All-Sky Monitor (ASM) count rate for 
the day of each Chandra observation (see Tab. 0J. Compar- 
ing the source flux against the ASM count rate, they are seen 
to be correlated, with no large offsets (see Fig. |3). 

To look for variability of the absorption lines between each 
epoch, in Fig. |4]we plot the measured equivalent width (or 
upper limit, if the line is not detected) for each observation vs 
the 0.5-10 keV unabsorbed source flux (see Table|4]l. The er- 
rorbars indicate la uncertainties in the equivalent width. The 
dashed lines indicates the weighted mean (weighted by the 
uncertainties) of the equivalent widths of the detected lines. 
Generally, the equivalent widths of the lines are consistent 
between the observations. However, we note that the O VIII 
line does seem to display variability, although the variability 
is not highly significant and does not appear to be correlated 
with the source flux (see Fig.|4]i. The observed wavelength of 
the O VIII line is also found to be consistent with the theoret- 
ical wavelength and not blueshifted. 

We calculate the x 2 value for a fit of the equivalent widths 
for the O VIII line to their weighted mean to assess the sig- 
nificance of this variability. For the non-detection in Obsld 
1021, we use the equivalent width determined from fitting a 
Gaussian fixed at the wavelength of the line in the calcula- 
tion. This equivalent width is found to 2.9 ±2.9 mA. We get 
X 2 = 14.1 (3 degrees of freedom) when fitting the equivalent 
widths to their weighted mean. For the hypothesis that the 
equivalent widths are constant, this corresponds to a proba- 
bility for achieving a higher x 2 value of 0.003. However, we 
note that we have searched 5 lines to find variability of this 
level, and thus this increases the probability to 0.015. This 
is equivalent to the line being variable at the 2.4<r level of 
confidence. Thus, while it is not a highly significant result, 
low-level variability cannot be ruled out for this line. 

In comparing equivalent width measurements from multi- 
ple epochs, it is important to consider whether any systematic 
effects can cause variability in the lines. Importantly, neither 
the first two observations (Obsld 1021 and 1022) were oper- 
ated in continuous-clocking mode, and thus both are affected 
by pileup. It is not clear how pileup alters equivalent width 
measurements, though we note that as the O VII and O VIII 
lines are close together in wavelength they should be affected 
similarly. Any effect must be fairly small, as the O VIII line is 
seen to vary whereas the O VII is not. 

5. DISCUSSION 

We have analyzed four Chandra HETGS observations of 
the neutron star low-mass X-ray binary 4U 1820-30 and de- 
tect a variety of absorption lines in the spectra. The first two 
archival observations suffer from pile-up, and thus measure- 
ments of lines from these observations may not be robust. 
However, the two new observations presented here do not suf- 
fer from this problem and have an increased sensitivity com- 
pared to the previous observations, allowing for a significantly 
improved study of the absorption lines in this system. 

To investigate the nature of these lines, we compare the line 
equivalent widths between the observations. If the line is un- 
saturated and on the linear part of the curve of growth, then 
the equivalent width is a direct measure of the absorbing col- 
umn. If the equivalent width does not vary between observa- 
tions, this would therefore indicate that the absorbing column 
has remained unchanged between the observations, as would 
be expected for absorption lines associated with the ISM. The 
majority of the absorption lines present in 4U 1820-30 are 
consistent with remaining constant between the observations, 



suggesting they are interstellar in origin. However, if the lines 
are saturated, then even if the column density changes signif- 
icantly (for instance in response to changes in the source) the 
equivalent width would remain more or less constant. Simi- 
larly, in the saturation regime even if the ionizing flux changes 
by a factor of 2 (as observed between Obsld 6633 or 6634) the 
equivalent width and ion relative abundance fractions may not 
change by a similar amount. However, even though we do not 
observe the a — (3 — 7 sequence for each line, upper limits on 
(3/ a indicate that the lines are not highly saturated. 

We find that there may be a low level of variability (2.4-ct 
significance) in the O VIII Lya line, which may indicate it is 
associated with the source. If the absorption is local to the 
source, this variability could be, for example, due to changes 
in mass accretion rate leading to changes in mass outflow. If 
one naively assumes that the continuum flux traces the mass 
accretion rate one would therefore see correlated variability 
between the source flux and the line equivalent width. Alter- 
natively, if the properties of the absorbing column are approx- 
imately constant and the ionizing flux changes, the fraction of 
oxygen in O VIII should change, which could lead to corre- 
lated variability. However, there appears to be no such cor- 
relation observed (see Fig. @|i though more observations are 
required to be conclusive. 

Additionally, if the gas in the system responds to changes 
in the ionizing flux more slowly than the inherent changes in 
the ionizing flux (i.e. if the recombination time is long), this 
would have the effect of greatly dampening the amplitude of 
variations seen in the absorption relative to the incident flux. 
For O VIII originating within the system, the size of the ab- 
sorbing region is < 10 n cm. The column density we observe 
to be about 10 2<) cm" 2 , thus the density is n e ^10 9 cm -3 . The 
recombination rate coefficient for O VIII is a = 1.5 x 10~ 12 
cm 3 s _I at 10 6 K, so the recombination timescale = l/(n e a) 
^100 s. This is significantly shorter than the timescale be- 
tween the observations and thus should not be a factor. 

If the absorption lines were associated with a disk wind, 
then blueshifts would be expected. We find that the wave- 
length of the absorption lines are consistent with their rest 
wavelengths (after accounting for the absolute instrumental 
calibration uncertainty of 0.05%), as expected for interstellar 
fines. However, we note that in the dipping X-ray binaries, 
absorption lines that are local to the source are not seen to be 
blueshifted and are thought to be associated with the accre- 
tion disk corona rather than an accretion disk wind. Similarly, 
here the lack of blueshifts could ju st be explained by abso rp- 
tion from an accretion disk corona (Diaz Tri go et aLl l2006). 

To test both possible origins for the lines, we first calculated 
the equivalent hydrogen column densities (implied by the ob- 
served equivalent widths) for a variety of Doppler parameters, 
and find that they are consistent with the interstellar origin for 
the lines. Additionally we perform photoionization model- 
ing of absorption by material local to the system. We find 
that photoionization can reproduce the observed spectra well. 
However, to do this requires a very low filling factor, but this 
could possibly be explained by absorption from dense blobs. 
We conclude that while the Ne II and Ne III lines are produced 
in the ISM, the origin of the Ne IX, O VII and O VIII lines re- 
mains unclear, with both the interstellar and local absorption 
remaining possible. The key to determining the origin of these 
lines may lie in variability, thus additional sensitive observa- 
tions of 4U 1820-30 are needed to confirm variability in the 
O VIII line, and search for low-level variability in other lines. 
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TABLE 1 

Chandra OBSERVATIONS OF 4U 1820-30 



ObsID Observation Date Grating/Detector Exposure (ks) 



98 2000 Mar 10 LETG/HRC 15.12 

1021 2001 Jul 21 HETG/ACIS 9.70 

1022 2001 Sep 12 HETG/ACIS 10.89 

6633 2006 Aug 12 HETG/ACIS 25.18 

6634 2006 Oct 20 HETG/ACIS 25.13 



NOTE. — ObsID 6633 and 6634 are the new observations presented in this work. 



TABLE 2 

Absorption line properties in 4U 1 820-30 



ObsID 


Ion 


Transition 


Theoretical 
(A) 


Observed 
(A) 


Equivalent width 
(mA) 


FWHM 

(mA) 


(io- 4 


Flux 
photons cm -2 s _1 ) 


1021 


NelX 
NeHI 
Ne II 
O VII 
O VIII 


ls 2 -ls2p 
ls-2p 
ls-2p 

Is 2 — ls3p 
ls-2p 


13.447 
14.508 
14.608 
18.629 
18.967 


13.455 ±0.003 


7.3 ±1.8 

< 3.7 

< 8.1 

< 10.2 

< 10.8 


3.8!" 




4.3 ±1.0 

< 2.1 

< 4.7 

< 4.8 

< 4.8 


1022 


NelX 
NeHI 
Ne II 
O VII 
O VIII 


ls 2 -ls2p 
ls-2p 
ls-2p 

ls 2 -ls3p 
ls-2p 


13.447 
14.508 
14.608 
18.629 
18.967 


13.448 ±0.005 
18.966 ±0.007 


6.4±2.2 

< 3.5 

< 3.6 
<6.6 

22.7 ±6.3 


19.3 ±10.9 
33.4±14.1 




5.2 ±1.8 

< 2.7 

< 2.8 

< 6.0 
1.4 ±0.4 


6633 


NelX 
NeHI 
Ne II 
O VII 
O VIII 


ls 2 -ls2p 
ls-2p 
ls-2p 

ls 2 -ls3p 
ls-2p 


13.447 
14.508 
14.608 
18.629 
18.967 


13.450 ±0.002 
14.505 ±0.006 
14.611 ±0.002 
18.634 ±0.006 
18.972 ±0.006 


4.0 ±0.9 
3.5 ± 1.3 

5.1 ±1.2 
9.8 ±3.0 
6.4 ±2.9 


1 4 +g ' ; 
18.6±13.5 

u ' u '-0.07 

24.0 ±11.0 

7 1 +2L2 




2.9 ±0.7 
2.5 ±0.9 
3.5 ±0.8 
6.3 ±2.0 
4.1 ± 1.8 


6634 


NelX 
NeHI 
Ne II 
O VII 
O VIII 


ls 2 -ls2p 
ls-2p 
ls-2p 

ls 2 -ls3p 
ls-2p 


13.447 
14.508 
14.608 
18.629 
18.967 


13.445 ±0.002 
14.508 ±0.003 
14.602 ±0.005 

1S.623_ 0001 
18.965 ±0.006 


5.0 ±0.8 

3.1 ±1.1 
5.4±1.2 
7.0±2.5 
18.4±4.5 


^ -1.4 
o 9+13.0 
-'■ z -3.2 

27. 8± 12.2 

o n? +I4 ' 4 

u - -0.02 

45.9 ±12.6 




4.9 ±0.8 

3.0 ±1.0 

5.1 ± 1.2 
6.1 ±2.3 
16.0 ±4.0 


Note. 


— See text for details on spectral fitting procedure. 


Theoretical line wavelengths for Ne IX, O VII, and O VIII are 


from Verner. Verner, & Ferland 



1 1996). The line wavelengths for Ne II and Ne III are the best observationally determined values from Juett et al. 12006). Uncertainties and upper limits are 
lcr confidence limits. Where the line was not detected and we determine upper limits, the wavelength was fixed at the theoretical value, and the width of the 
Gaussian fixed below the instrument resolution. 



TABLE 3 

Line column densities for various Doppler parameters 



ObsID 


Line 


Linear curve of growth 


b = 


= 50 km s _1 


b = 


100 km s" 1 


b = 


= 200 km s"" 1 






N io „ 


N H 


N ion 


N H 


N itm 


N H 




N H 






(10 16 cm 4 ) 


(10 20 cm" 2 ) 


(10 16 cm" 


" 2 ) (10 20 cnT 2 ) 


(10 16 cm" 


2 ) (10 20 cnT 2 ) 


(10 16 cm" 


2 ) (10 20 cm" 2 ) 


6633 


NelX (13.447 A) 


0.3 


0.4 


0.6 


0.7 


0.4 


0.5 


0.4 


0.4 




Ne III(14.508A) 


1.8 


7.1 


2.6 


12.5 


2.1 


8.9 


1.9 


7.8 




NeII(14.608A) 


4.4 


7.1 


8.3 


12.6 


5.6 


8.9 


4.9 


7.8 




OVII(18.629A) 


2.2 


0.5 


9.7 


2.1 


3.3 


0.8 


2.6 


0.6 




O VIII (18.967 A) 


0.5 


0.6 


0.7 


2.1 


0.6 


0.9 


0.5 


0.7 


6634 


NeIX(13.447A) 


0.4 


0.5 


0.9 


1.1 


0.6 


0.7 


0.5 


0.6 




Ne III (14.508A) 


1.6 


7.1 


2.2 


13.4 


1.8 


9.0 


1.7 


7.9 




NeII(14.608A) 


4.6 


7.1 


9.4 


13.4 


6.0 


9.1 


5.2 


8.0 




OVII(18.629A) 


1.6 


0.6 


3.2 


1.3 


2.1 


0.8 


1.8 


0.7 




O VIII (18.967 A) 


1.4 


0.8 


6.4 


3.2 


2.9 


1.5 


1.9 


1.0 



NOTE. — See text for details on determining column densities. We only determine column densities for ObsIDs 6633 and 6634. The equivalent hydrogen 
column density (Nh) is determined using the ISM abundances from Wilms et al. 1 2000). 
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TABLE 4 
Spectral fits to 4U 1820-30 



ObsID 




Blackbody 


Power-law 




0.5 — 10 keV unabsorbed flux 


RXTEI ASM count rate 




kT (keV) 


Normalization 


Spectral index 


Normalization 




(10 -9 erg cirT 2 s~') 


(counts s ) 


1021 


0.72±0.07 


(3.5 ±0.6) x 10~ 3 


1.67 ±0.02 


0.85 ±0.02 


0.87 


5.8±0.2 


14.0±0.3 


1022 


9fi+ 015 
-0. 1 1 


(4.5 ± 1.0) x 10~ 3 


1.46 ±0.02 


1.13 ±0.02 


0.88 


9.7 ±0.3 


24.3 ±0.3 


6633 


1.16±0.02 


(2.1 ±0.1) x 10~ 2 


1.76 ±0.03 


1.12±0.01 


0.94 


8.3 ±0.2 


21.3±0.6 


6634 


1.20 ±0.02 


(3.2 ±0.2) x Kr 2 


1.74 ±0.02 


1.50±0.01 


1.24 


11.7±0.3 


27.7±1.0 



NOTE. — In all fits the column density was fixed at the cluster value of Nh = 1.5 x 10 21 cm -2 . We used the XSPEC blackbody model named bbody, where the 
normalization is defined as L39/(D10) 2 , where L39 is the luminosity in units of 10 3 ' erg s and D10 the distance to the source in units of 10 kpc. The power-law 
normalization is defined at 1 keV, with units of photons keV -1 cm -2 s . The spectral fits are to the HEG 1st order spectrum (see text for details). For comparison 
with the source flux, the RXTEI ASM count rate for 4U 1 820—30 is also given. 




Wavelength [Angstrom] 



Wavelength [Angstrom] 




_ i i i i i i 



13.5 14 14.5 

Wavelength [Angstrom] 



18.5 19 19.5 

Wavelength [Angstrom] 



FIG. 1 . — 13 - 15A and 18 — 20A segments of the combined first-order MEG spectra of 4U 1820—30. For illustrative purposes the data have been rebinned by 
a factor of 2. The top panels show these regions in Obsld 6633, and the bottom panels show these regions in Obsld 6634. The data is in black, and the \a error 
bars are in blue. The best-fitting model is in red. The continuum is fit by a power-law modified by neutral ISM absorption edges, where appropriate. Absorption 
lines are modeled by a Gaussian. Table[2]gives the line properties. 
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FIG. 3. — 0.5-10 keV unabsorbed source flux determined from fitting to the HEG spectra is shown versus the 1-day average RXTE All-Sky Monitor count rate 
on the day of each observation. 



10 



Cackett et al. 



«; o 



O 

CO 



o 

C\2 



CO 
CD 

CD 



CV2 



m - 



\ 


/i . i 


. \ 


i 1 i . 

VII : 

' i i 


: \ 


i 1 i 1 


H 1 1 — = 

VIII : 

r + 1 

1,1" 


\ 


— +- - - 

'1 , I , v 


H 1 1 

Ne II ~ 

----+---■ 

' i , i 


\ 


1 

'1 , 1 


i 


H 1 1 

Ne III " 

—+--: 

' i , i 


1 1 1 1 

; + 

+ ._ _| 


H 1 1 

Ne IX - 

7 1 : 



0.6 0.8 



1 



Flux (10 8 erg cm a s 



1.2 
- 1 ) 



FIG. 4. — Equivalent widths (mA) of the various absorption lines in 4U 1820-30 versus 0.5-10 keV unabsorbed source flux in units of 10~ 8 erg cirT 2 
. Dashed lines indicate the weighted mean of the measured equivalent widths. Arrows indicate 95% confidence upper limits and errorbars indicate 1 
uncertainties. 



